Intermittent Switching between Potential Flow and Turbulence in Superfluid Helium 

at mK Temperatures 
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Superfluid flow around an oscillating microsphere is investigated at temperatures down to 25 mK. 
Stable laminar flow below a critical velocity and turbulence at large drives are found to be separated 
below 0.5 K by an intermediate range of driving forces where the flow is unstable, intermittently 
switching between laminar and turbulent phases. We have recorded time series of this switching 
phenomenon and have made a statistical analysis of the switching probability. The mean lifetime of 
the turbulent phases grows with increasing drive and becomes infinite at a critical value. Stability 
of the laminar phases above the critical velocity is limited by natural background radioactivity or 
cosmic rays. 

PACS numbers: 67.40.Vs, 47.27.Cn 



Turbulence in a classical liquid is based on the nonlin- 
earity of the Navier-Stokes equation of motion H while 
turbulence in superfluid helium, which also has been a 
very active area of research for several decades, is de- 
scribed by two-fluid hydrodynamics of the normal and 
superfluid components of the liquid ppl] . Because of the 
quantum mechanical properties of the liquid, turbulence 
in superfluid helium consists of a tangle of individual vor- 
tex lines having the same quantized circulation h/m^, 
where h is Planck's constant and nn is the mass of a 
helium-4 atom. In this respect, superfluid turbulence is 
different from classical turbulence but it has also striking 
similarities like a Kolmogorov energy spectrum in the in- 
ertial range g. Both cases represent different branches 
of the general phenomenon of fluid turbulence. A com- 
pletely new aspect arises when superfluid turbulence is 
considered in the limit of very low temperatures [|J . In 
this case all what is left of the normal fluid component 
is a very dilute gas of ballistically propagating phonons 
and therefore turbulence in the pure superfluid can be 
investigated. While a number of experiments have been 
performed which demonstrate the production of individ- 
ual vortex lines both in superfluid He and in superfluid 
3 He |j) and the decay of the vortex tangle J7|, the tran- 
sition to turbulence in the flow of the superfluid around 
a macroscopic body so far has not been investigated at 
these very low temperatures. 

In our present work, we have studied superfluid flow 
around an oscillating microsphere at temperatures down 
to 25 mK. Laminar and turbulent flow can be identified 
by the drag force on the sphere, which in case of laminar 
flow is linear and determined by ballistic phonon scatter- 
ing, while for turbulent flow, which exists above a criti- 
cal velocity, the drag force is nonlinear and much larger 
[pj. Here, we report on the observation of an interme- 
diate regime located between stable laminar flow below 
the critical velocity and stable turbulence above a criti- 
cal driving force, where the flow is found to be unstable, 



switching intermittently between laminar and turbulent 
phases. We have measured time series of this switch- 
ing phenomenon at different constant temperatures and 
driving forces. We find, firstly, the mean lifetime of the 
turbulent phases to grow with increasing drive and to 
diverge at a critical value. Secondly, the stability of the 
laminar phases above the critical velocity is analyzed and 
the probability for switching to turbulence is obtained. 
Surprisingly, we find the long term stability of the lam- 
inar phases close to the critical velocity to be limited 
ultimately by natural radioactivity or cosmic rays. 

Our experimental technique has been described in ear- 
lier work S and is only briefly summarized here as 
follows. The microsphere (radius R — 0.1mm, mass 
m = 27 /ig) is made of ferromagnetic SmCos. Therefore, 
we can use superconducting levitation, i.e., the repulsion 
of a permanent magnet from a superconducting surface, 
to let it float in the middle between the superconducting 
electrodes of a parallel plate capacitor made of niobium 
(spacing d — 1 mm, diameter 4 mm) without any me- 
chanical suspension. Lateral stability is provided by flux 
lines trapped in the superconductor. The sphere carries 
an electric charge q of about 1 pC when it levitates, be- 
cause we apply a voltage of several hundred volts to the 
bottom electrode when cooling the capacitor through the 
critical temperature T c = 9.2 K of niobium. Vertical os- 
cillations of the sphere about its equilibrium position can 
now be excited by applying a small ac voltage U ac ranging 
from 0.1 mV up to 20 V at the resonance frequency of the 
sphere ~ 120 Hz. These oscillations induce an ac current 
qv/d in the electrodes, where v is the velocity amplitude, 
which is detected by an electrometer. The damping of 
the oscillations is very small, Q- values in excess of one 
million are achieved if the capacitor is evacuated M . We 
then fill the capacitor with pure superfluid helium |1C| ] 
and measure the velocity amplitude at resonance as a 
function of the driving force F = qU ac /d at constant 
temperature between 25 mK and 1 K. 



As reported earlier S there is a linear regime v(F) at 
small amplitudes, see Fig. 1. This indicates a linear drag 
force X(T)v where the drag coefficient A oc T 4 is deter- 
mined by ballistic scattering of thermal phonons which 
can be considered as a dilute gas of excitations whose 
density rapidly decreases as T 4 . The flow of the super- 
fluid around the sphere is frictionless potential flow. At 
large driving forces, however, the sphere creates vortic- 
ity in the superfluid. Now the drag force is large and 
nonlinear, see Fig. 1. We find the drag force to be given 
by Fd = jv 2 — F , where 7 = C£>pnR 2 /2 is the classical 
turbulent drag coefficient (p is the density of the liquid 
and cd ~ 0.4 for a sphere) and Fq is a constant. In 
contrast to a classical liquid where Fq — in our case 
turbulent drag (Fd > 0) exists only above a critical ve- 
locity v c — 19.4 mm/s. A similar drag force has recently 
been calculated for a cylinder in a two-dimensional dilute 
Bose-Einstein condensate, although the critical velocity 
is very much different in that case [O. 
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FIG. 1. Velocity amplitude of the oscillating sphere as a 
function of the driving force at 300 niK. At small drives the 
flow is laminar and the linear behavior is determined by ballis- 
tic phonon scattering (•). The shaded area marks the regime 
of instability above the critical velocity v c . At larger driving 
forces we observe stable nonlinear turbulent drag (■), where 
the solid line is a fit of the quadratic drag force Fd, see text. 

Most interesting is our new observation of an interme- 
diate regime where neither potential flow nor turbulence 
are stable but instead the flow switches between both 
patterns intermittently [O] . This regime of instability is 
observed only below 0.5 K and extends from the critical 
velocity v c up to a critical driving force where enough 
power is being delivered to the superfluid turbulence to 
become stable. We have recorded time series (up to 36 
hours per run) of this switching phenomenon at various 
driving forces and temperatures. An example is shown in 
Fig. 2 where a section of 17 minutes is shown for three dif- 
ferent driving forces at 300 mK. The amplitude switches 
between a low level corresponding to turbulent drag and 
an exponential recovery of the laminar level correspond- 



ing to phonon scattering. If the laminar phase lasts long 
enough the stationary value for laminar flow is reached 
(see upper time series). When the laminar flow breaks 
down a rapid drop to the turbulent level occurs. It is 
obvious from Fig. 2 that with increasing driving force the 
lifetimes of the laminar phases shrink while those of the 
turbulent phases grow. 
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FIG. 2. Intermittent switching of the velocity amplitude at 
300 mK. The low level corresponds to turbulent flow while the 
increase occurs during a laminar phase. The driving forces 
(in pN) are 47, 55, 75 (from top to bottom). 
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FIG. 3. Statistical analysis of the lifetimes t of the turbu- 
lent phases, a) The exponential distribution P(t) at a given 
temperature and driving force, b) The slopes v measured 
at different temperatures as a function of the turbulent drag 
force. The solid line is a fit of a fourth power divergence of 
l/u at a critical force of 54.6 pN. 

In the following we discuss the results of a statistical 
analysis of these time series based on reliability theory 
]13[ . A typical distribution of the lifetimes of the tur- 
bulent phases is shown in Fig. 3a where the number P 
of turbulent phases having a lifetime longer than a given 
threshold t is plotted on a logarithmic scale versus t. We 
always find an exponential distribution P(t) oc exp(— i/t) 
with 1/v being the mean lifetime. This implies a prob- 
ability for the breakdown of the turbulent phase ("fail- 
ure rate") —dhxP/dt = v being constant in time (i.e 



no history dependence). The values of v are indepen- 
dent of temperature and collapse to an universal drive 
dependence if the strongly temperature dependent lam- 
inar drag Xv (where v = vt is the velocity amplitude 
of the turbulent phase) is subtracted from the external 
driving force. Obviously, v depends only on the turbu- 
lent drag force F — Xv t and vanishes at a critical drive 
approximately as the fourth power of the distance from a 
critical value, see Fig. 3b. At the critical drive the power 
delivered to the turbulent supcrfluid is 0.6 pW which cor- 
responds approximately to the production of one vortex 
ring of the size of the sphere per half period of the oscil- 
lation. 
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FIG. 4. Statistical analysis of the laminar phases, a) The 
Gaussian distribution P(Av) oc exp( — (Av/vw) 2 ) of the ve- 
locity increase Av at a given temperature and driving force, 
b) The fitting parameter vw is shown to be independent of 
temperature and driving force. 

Stability of the laminar phases can be analyzed in two 
ways, either the distribution of the signal height Av (ve- 
locity at breakdown minus turbulent velocity vt) is con- 
sidered or the distribution of the lifetimes. Both vari- 
ables are related by the exponential recovery Av(t) = 
Aw maa; (l — exp(— t/r)), where Av max is the difference 
between the stationary laminar velocity amplitude F/X 
and vt, and r = 2m/ X is the time constant. In Fig. 4a 
the distribution of the signal height Av is shown to be 
of the Gaussian form P(Av) oc exp(— (Av/vw) 2 ) where 
vw is a fitting parameter. This implies a failure rate (per 
unit velocity increment) —d\nP/dAv = 2Avjv w to be 
proportional to the signal height Av. The corresponding 
probability density of Av (—dP/dAv) is called a Weibull 
distribution. The fitting parameter vw ~ 4.8 mm/s is 
found to be completely independent of both tempera- 
ture and driving force, see Fig. 4b. We expect it to de- 
pend on the size or the surface roughness of the sphere 
but experiments with different spheres remain to be per- 
formed in the future. Because Av(t) saturates for Or, 
it is more informative in this case to analyze the life- 
time distribution, in particular near v c where the lam- 
inar phase becomes stable. From P(t) = P(Av(t)) oc 
exp(— (Av(t)/vw) 2 ) we conclude that for i ^ t P{t) 
must be constant and therefore the failure rate (per unit 



time) is -dhiP/dt = v w 2 d(Av) 2 /dt = 0. This simply 
means that if the laminar phase has survived for a time 
t 3> r it will be stable indefinitely although the velocity 
is clearly above v c . Our experimental results, however, 
demonstrate that the lifetime is finite: \nP(t) slowly de- 
creases at large t and from a straight line fit to the data 
in Fig. 5 we obtain a mean lifetime of 25 minutes. We 
can rule out mechanical vibrations of the cryostat as the 
source of this instability: jumping on the floor, slamming 
the door, even transferring liquid helium to the cryostat 
had no effect. Therefore, another mechanism must be 
responsible. In fact, placing a small radioactive source 
60 Co (74kBq) near the cryostat has a dramatic effect: 
the mean lifetime of the laminar phase is reduced to 3.0 
minutes. We have measured the dose rate of the source at 
the position of the measuring cell inside the cryostat (tak- 
ing into account a measured 20% loss in the dewar walls) 
to be 440nGy/h (±5%). Comparing this value with a 
measured dose rate due to natural background radiation 
in our laboratory of 50nGy/h (± 10%), which is typical 
for our area, we obtain an increase of the dose rate due 
to the source by a factor of (440 + 50)/50 = 9.8. Within 
the error bars this compares well with the measured re- 
duction of the mean lifetime of the laminar phases by a 
factor of 25/3.0 = 8.3. This effect may be attributed to 
the creation of ions in the superfluid, which produce lo- 
cal vorticity either during the creation and recombination 
processes or when the ions are accelerated by the electric 
field which exists in our measuring cell |L4|. This vor- 
ticity then initiates the breakdown of the potential flow 
if the velocity amplitude of the sphere is larger than the 
critical velocity v c . We therefore conclude that the life- 
time of the metastable laminar phases which we observe 
above v c is only limited by natural background radioac- 
tivity. 
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FIG. 5. Distribution of the lifetimes of the laminar phases 
slightly above the critical velocity v c . For t 2> t = 31 s (at 
300 mK) the mean lifetime of 25 minutes is reduced by the 
radioactive source to 3 minutes. 

In summary, we have found that the transition from 



potential flow to turbulence around a sphere in super- 
fluid helium below 0.5 K occurs by intermittent switch- 
ing between both flow patterns. We understand now 
why above 0.5 K the intermittent switching changes into 
the hysteretic behavior observed earlier m: the velocity 
increases Av become very small at higher phonon drag 
which implies a very low failure rate of the laminar phase. 
But if it fails (i.e. when the driving force F is largely in- 
creased) the following turbulent phase is stable because 
the critical drive is exceeded. 

At present we have no theoretical model which could 
describe our experimental results for the failure rates of 
the laminar and turbulent phases. But because the tran- 
sition to turbulence in superfluid helium at mK tempera- 
tures is a new phenomenon, we hope that our results will 
stimulate further theoretical progress. 
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